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Newtonian versus relativistic nonlinear cosmology 



Hyerim Noh* and 


Abstract 


Both for the background world model and its linear perturbations New¬ 
tonian cosmology coincides with the zero-pressure limits of relativistic cos¬ 
mology. However, such successes in Newtonian cosmology are not purely 
based on Newton’s gravity, but are rather guided ones by previously known 
results in Einstein’s theory. The action-at-a-distance nature of Newton’s 
gravity requires further verification from Einstein’s theory for its use in the 
large-scale nonlinear regimes. We study the domain of validity of the New¬ 
tonian cosmology by investigating weakly nonlinear regimes in relativistic 


cosmology assuming a zero-pressure and irrotational fiuid. We show that, 
first, if we ignore the coupling with gravitational waves the Newtonian cos¬ 
mology is exactly valid even to the second order in perturbation. Second, 
the pure relativistic correction terms start appearing from the third order. 
Third, the correction terms are independent of the horizon scale and are 
quite small in the large-scale near the horizon. These conclusions are based 
on our special (and proper) choice of variables and gauge conditions. In a 
complementary situation where the system is weakly relativistic but fully 
nonlinear (thus, far inside the horizon) we can employ the post-Newtonian 
approximation. We also show that in the large-scale structures the post- 
Newtonian effects are quite small. As a consequence, now we can rely on 
the Newtonian gravity in analyzing the evolution of nonlinear large-scale 
structures even near the horizon volume. 

1. Introduction: In order to interpret results from Einstein’s gravity theory properly 
we often need corresponding results in Newton’s theory. On the other hand, in order to use 
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results from Newton’s gravity theory reliably we need confirmation from Einstein’s theory. 
The observed large-scale structures show nonlinear processes are working. Currently, studies 
of such structures are mainly based on Newtonian physics in both analytical and numeri¬ 
cal approaches. One may admit its incompleteness as the simulation scale becomes large 
because, first, Newton’s gravity is an act ion-at-a-distance, i.e., the gravitational influence 
propagates instantaneously thus violating causality. Second, Newton’s theory is ignorant of 
the presence of horizon where the relativistic effects are supposed to dominate. One other 
reason we may add is that Einstein’s gravity apparently has quite different structure from 
Newton’s one. The causality of gravitational interactions and consequent presence of the 
horizon in cosmology are naturally taken into account in the relativistic gravity theory. 

In the literature, however, independently of such possible shortcomings of Newton’s 
gravity in the cosmological situation, the physical size of Newtonian simulation, in fact, has 
already reached the Hubble horizon scale. Common excuses often made by people working 
in this active field of large-scale numerical simulation are, hrst, in the small scale one may 
rely on Newton’s theory and, second, as the scale becomes large the large-scale distribution 
of galaxies looks homogeneous, //the deviation from homogeneity is small (linear) Einstein’s 
gravity gives the same result as the Newtonian one. Presence of large-scale homogeneity, 
although difficult to verify observationally, is in fact a crucially important assumption in 
currently popular cosmology. In order to have proper conhrmation, however, we still need to 
investigate Einstein’s case in the nonlinear or weakly nonlinear situations. While the general 
relativistic cosmological simulation is currently not available, in this work, we will shed light 
on the situation by a perturbative study of the nonlinear regimes assuming zero-pressure 
and irrotational fluid in Einstein’s gravity. This allows us to investigate the similarity and 
difference between the two gravity theories in the weakly nonlinear regimes in cosmological 
situation. We will show that even to the second order in perturbations, except for the cou¬ 
pling with gravitational waves, Einstein’s gravity gives the same results known in Newton’s 
theory and the pure relativistic corrections appearing in the third order perturbations are 
independent of the horizon and are small. We also present a complementary approach using 
the post-Newtonian approximation which can handle weakly relativistic (thus, far inside 
the horizon) but fully nonlinear situation. We show that the hrst-order post-Newtonian 
corrections are again quite small. Thus, now our relativistic analysis assures that Newton’s 
gravity is practically reliable even in the weakly nonlinear regimes in cosmology. We set 
c = 1. 


2. Nonlinear equations: We start from the completely nonlinear and covariant equa¬ 
tions [1]. We need the energy conservation equation and the Raychaudhury equation. In a 
zero-pressure medium without rotation we have [1] 

jl fl6 = 0, 6 —6^ d'^^^d'ab + 4:TTGfl — A = 0, (1) 

o 

where A is the cosmological constant; jl is the energy density, 6 = 12“.^ is the expansion 
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scalar with Ua the fluid four-vector, and aat is the shear tensor; jl = is the covariant 
derivatives along u°‘. Tildes indicate the covariant quantities. By combining these equations 
we have 



— a^'^aob — AttG/x -|- A = 0. 


( 2 ) 


These equations are fully nonlinear and covariant. To the second and higher order pertur¬ 
bations we also need the momentum constraint part of Einstein’s equation. 

As the metric we take 


ds'^ = —a^ (1 -I- 2a) drf — 2a^ (3^adrjdx°‘ -|- (1 -|- 2ip) + 2')^a\(i dx'^dx^ , 


( 3 ) 


where a, 13, j and (p are spacetime dependent perturbed-order variables. Spatial indices 

/o\ 

of perturbed order variables are based on and a vertical bar indicates the covariant 
derivative based on We ignored the transverse vector-type perturbation and transverse- 
tracefree tensor-type perturbation variables. In this perturbation study we will consider 
the scalar-type perturbations up to third order in the flat Friedmann background without 
pressure. The vector-type perturbation has only a decaying solution in expanding medium. 
The presence of gravitational waves will cause couplings with the scalar-type perturbation 
to the second and higher orders in perturbations, see [2, 3]. The presence of gravitational 
waves can be regarded as a pure relativistic effect even to the linear order. 


3. Background world model: To the background order, we have fi = jj, and 9 = 3^ 
where a{t) is the scale factor, and an overdot indicates the time derivative based on the 
background proper-time t. Equation (1) gives 

fi 3-/i = 0, 3- -|- AnGjJi — A = 0. (4) 

a a 

This was first derived based on Einstein’s gravity by Friedmann in 1922 [4], and the New¬ 
tonian study followed later by Milne and McCrea in 1934 [5]. In the Newtonian context fi 
can be identified with the mass density g. 


4- Linear-order perturbations: To the linear-order perturbations in the metric and 
energy-momentum variables, we introduce 

f = fi -\- 6fi, 9 = 3 —h 69. (5) 

a 

To the linear order we identify 

6fi = 6q, 69 = -V • u, (6) 

a 

where 6g and u are the perturbed mass density and the peculiar velocity in Newtonian 
context. In all our relativistic (nonlinear) perturbation analyses we take the temporal 
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comoving gauge (which together with the irrotational condition gives «« = 0 for the fluid 
four-vector Ua) and the spatial 7 = 0 gauge [6]. As these gauge conditions fix the gauge 
modes completely, all the remaining variables are equivalently gauge-invariant to all orders 
in perturbations; for technical details, see Section VI.C of [6] where an explicitly gauge- 
invariant combination which is the same as the Sfi in the temporal comoving gauge and 
the spatial 7 = 0 gauge can be found in eq. (282). These gauge conditions and our choice 
of the perturbation variables are crucially important to make our conclusions. 

To the linear order the perturbed part of eq. (2) gives 

S + 2-6- ATTGfid = 0, (7) 

a 

which is the well known density perturbation equation in both relativistic and Newtonian 
contexts; we set S = 5n/fi. This equation was first derived based on Einstein’s gravity by 
Lifshitz in 1946 [7], and the Newtonian study followed later by Bonnor in 1957 [8]. 

It is curious to notice that in both the expanding world model and its linear structures 
the first studies were made in the context of Einstein’s gravity (Friedmann 1922; Lifshitz 
1946), and the much simpler and, in hindsight, more intuitive Newtonian studies followed 
later (Milne and McCrea 1934; Bonnor 1957). Perhaps these historical developments reflect 
that people did not have confidence in using Newton’s gravity in cosmology before the result 
was already known in, and the method was ushered by, Einstein’s gravity. It may be also 
true that only after having a Newtonian counterpart we could understand better what the 
often arcane relativistic analysis shows. It would be fair to point out, however, that the 
ordinarily known Newtonian cosmology (both for the background world model and its linear 
perturbations) is not purely based on Newton’s gravity, but is a guided one by Einstein’s 
theory [9]. In the cosmological context Newtonian gravity is known to be incomplete and 
inconsistent; these are due to lack of boundary condition at spatial infinity and the action- 
at-a-distance nature of Newton’s gravity. For the second-order perturbations, currently we 
only have the Newtonian result known in the literature. Thus, the result only known in 
Newton’s gravity still awaits confirmation from Einstein’s theory. Here, we are going to fill 
the gap by presenting the much needed relativistic confirmation to the second order and 
the pure general relativistic corrections start appearing from the third order [2, 3]. 

5. Second-order perturbations: Even to the second order we introduce perturbations 
as in eq. (5), and take the same identifications made in eq. (6). To the second order the 
perturbed part of eq. (2) gives [6, 2] 

S + 2-6 — 47rGja6 = —r— [aV • (5u)] H —;rV • (u ■ Vu), (8) 

a ot 

where the second-order terms are in the right hand side. Exactly the same equation also 
follows from Newton’s theory [10]. Although we identified the relativistic density and 
velocity perturbation variables we cannot identify a relativistic variable which corresponds 
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to the Newtonian gravitational potential to the second order [2]. This may not be surprising 
because Poisson’s equation indeed reveals the action-at-a-distance nature and the static 
nature of Newton’s gravity theory compared with Einstein’s gravity. In the Newtonian 
context eq. (8) is valid to fully nonlinear order. 


6. Third-order perturbations: Since the zero-pressure Newtonian system is exact to the 
second order in nonlinearity, all non-vanishing third and higher order perturbation terms 
in the relativistic analysis can be regarded as the pure relativistic corrections. We use the 
same identification made in eq. (6) to be valid even to the third order, and will take the 
consequent additional third order terms as the pure relativistic corrections. To the third 
order the perturbed part of eq. (2) gives [3] 


5 -I- 2-6 — inGfiS = —r— aV • (5u) H—-V • (u ■ Vu) 
a a^ ot a^ 

l|{a[2^u-V(A-‘X)].Vq-lv 


a? dt 
2 


(^u • V (V • u) + 4 [u • V (A-^x) 


3a^ 


(p u • Vu-uV • u 

V 3 

4" ■ VX - • u, 


?>a? 


where the last two lines are pure third-order terms with 

3 

X = 2LpV • u — u • -|- -A“^V • [u • V (V 99 ) -|- uA^?]. 


( 9 ) 


This extends eq. (8) to the third order. Notice that we need the behavior of (p to the 
linear order only; 99 is a perturbed part of three-space metric in eq. (3), related to the 
perturbed three-space curvature (in our comoving gauge), and dimensionless. The third- 
order correction terms in eq. (9) reveal that all of them are simply of v^-order higher than 
the second-order terms. Thus, the pure general relativistic effects are at least v^-order higher 
than the relativistic/Newtonian ones in the second order. To the linear order we have [11] 


0 = 0, 


( 10 ) 


thus 99 = C(x) with no decaying mode. To the linear order our p is related to the perturbed 
Newtonian potential (5<l> and the Newtonian peculiar velocity u as [2, 3] 

(p = —(5<h-|-dA“^V • u. (11) 


The temperature anisotropy of cosmic microwave background radiation gives [12, 13] 
6T 1 


:P 


10 


-5 


r 5— - 

Thus, 99 ~ 5 X 10“® in the large-scale limit near horizon scale. Therefore, to the third 
order, the pure relativistic corrections are independent of the horizon scale and depend on 
the strength of linear order curvature perturbation p only, and are small. 
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7. Discussion: In this work we show that Newtonian cosmological perturbation equa¬ 
tions remain valid in all cosmological scales including the super-horizon scale to the second 
order. We assumed a zero-pressure irrotational fluid and ignored the coupling with grav¬ 
itational waves. The pure general relativistic correction terms start appearing from the 
third order. The third order correction terms involve only (p which is independent of the 
horizon scale and is small in the large scale limit near horizon. Therefore, one can now 
use the large-scale Newtonian numerical simulation more reliably as the simulation scale 
approaches and even goes beyond the horizon. All our results include the cosmological 
constant thus relevant in currently favoured cosmology. 

The referee has raised a couple of interesting observations that our conclusions do not 
refer to the averaging procedure [14], and the pure relativistic corrections start appearing 
at third order do not depend on the physical scale and on the averaging procedure. Indeed, 
our relativistic-Newtonian correspondence and the pure relativistic correction terms do 
not depend on scales nor on averaging procedure. We have reached our conclusions by 
comparing the exact Newtonian equations with the relativistic ones perturbed to the second 
and third orders without taking any averaging procedure. Thus, our relativistic-Newtonian 
correspondence to the second order and pure relativistic correction terms to the third order 
are independent of the averaging procedure. Notice, however, that we have achieved our 
result by choosing special (and proper) variables in certain (spatial and temporal) gauge 
conditions where all the variables have corresponding unique gauge-invariant combination 
of variables. 

The independence of the third order pure relativistic correction terms from the scale 
(compared with the second-order terms) is a sure surprise of our result. However, we would 
like to point out that our pure relativistic correction terms in eq. (9), certainly depend on our 
identification of the relativistic gauge-invariant combination of variables as the Newtonian 
ones to the third order made in eq. (6); this point was emphasized above eq. (9). Thus, 
if we take other identification (of the relativistic variables and gauges) as the Newtonian 
ones we could end up with correction terms which differ from our result. Based on our 
successful and clear identification with exact relativistic-Newtonian correspondence to the 
second order we believe (therefore, propose) the same identification to be valid to the third 
order, and suggest the third order correction terms in eq. (9) as the pure general relativistic 
effects (based on our identification of the variables). 

The roles of tensor-type perturbation (gravitational waves) are studied in [2, 3]; vector- 
type perturbation (rotation) is not important because it always decays in the expanding 
phase. Why Newtonian cosmology, despite its action-at-a-distance nature, still gives the 
same relativistic results even to the second-order perturbation in all scales, leaves room for 
further clarification. Also, it would be interesting to find cosmological situations where the 
pure general relativistic correction terms in eq. (9) could have observationally distinguish¬ 
able consequences. 

Consistency of the Newtonian (nonlinear) cosmology with the Newtonian limit of the 
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post-Newtonian approximation of general relativity was also reported in [15]. In fact, it 
is well known that the Newtonian hydrodynamic equations naturally appear in the zeroth 
post-Newtonian order of Einstein’s gravity [16]. In [15] it was shown that it is essential 
to keep the magnetic part of Weyl tensor in order to properly recover even the Newtonian 
limit in the post-Newtonian approach. In making our proof of the relativistic-Newtonian 
correspondence to the second order we assumed irrotational and zero-pressure conditions 
but have not imposed any condition on the magnetic part of the Weyl tensor; for a study 
based on the covariant equations ignoring the latter quantity, see [17]. In fact, the magnetic 
part of Weyl tensor does not vanish even to the linear order in perturbations: this quantity 
valid to the second order is presented in eq. (96) of [6]. 

Our nonlinear perturbation approach is applicable to fully relativistic regimes includ¬ 
ing the super-horizon scales and the early universe. However, it is limited to the weakly 
nonlinear situations where the nonlinearity is supposed to be small. A complementary 
approach in handling the large-scale nonlinear evolution in Einstein’s gravity is the post- 
Newtonian approximation. The post-Newtonian approach assumes w/c-expansion with 
GM/{Rc^) ~ v'^/c^ -C 1. Whereas our perturbation approach is applicable in fully rel¬ 
ativistic regime assuming weak nonlinearity, the post-Newtonian approach is applicable in 
fully nonlinear regime assuming weak (relativistic) gravity and slow motion. Thus, whereas 
the perturbation approach is applicable in all scales assuming weak nonlinearity, the post- 
Newtonian approach is applicable to fully nonlinear stage but only inside the horizon. 
Therefore, these two approaches are complimentary in the research of large-scale cosmic 
structures. Recently, we have extended Chandrasekhar’s first-order post-Newtonian hy¬ 
drodynamic approximation [16] to cosmological situation [18]. In [18] we show that the 
hrst-order post-Newtonian correction terms are of order 


GM 


10 “® 


10 


-4 


(13) 


compared with the Newtonian terms. Thus, although there could appear secular effects 
due to time-delayed propagation of gravity, the relativistic corrections are quite negligibly 
small similarly as our third-order pure relativistic correction terms in the weakly nonlinear 
regime. 

Therefore, our weakly nonlinear perturbation study and the fully nonlinear post- 
Newtonian study assure that in the current stage of the large-scale structure the Newtonian 
hydrodynamic equations are quite sufficient and reliable in handling the dynamics. How¬ 
ever, since we have not identihed the relativistic variable which corresponds to Newtonian 
potential to the second order, the Newtonian equations are not supposed to be reliable 
where the gravitational potential has an important role, like the gravitational lensing. 
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